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ABSTRACT 

We present a study of the composition of gas and dust in the Large and Small Magellanic Clouds 
(LMC and SMC, together - the MCs) as measured by UV absorption spectroscopy. We have measured 
P H and Fe H along 85 sightlines toward the MCs using archival FUSE observations. For 16 of those 
sightlines, we have measured Si H, Cr H, and Zn H from new HST COS observations. We have 
combined these measurements with H I and H 2 column densities and reference stellar abundances from 
the literature to derive gas-phase abundances, depletions, and gas-to-dust ratios (GDRs). 80 of our 84 
P measurements and 13 of our 16 Zn measurements are depleted by more than 0.1 decades, suggesting 
that P and Zn abundances are not accurate metallicity indicators at and above the metallicity of the 
SMC. The maximum P and Zn depletions are the same in the MW, LMC, and SMC. Si, Cr, and Fe 
are systematically less depleted in the SMC than in the MW or LMC. The minimum Si depletion in 
the SMC is consistent with zero. Our depletion-derived GDRs broadly agree with GDRs from the 
literature. The GDR varies from location to location within a galaxy by a factor of up to 2 in the 


LMC and up to 5 in the SMC. This variation is 
diffuse neutral phase of the interstellar medium. 
Subject headings: Galaxies: abundances, galaxies: 

1. INTRODUCTION 

From different types of observations of dust in the 
Milky Way (MW) and other galaxies, we know that 
the gas-to-dust ratio (GDR) varies within and between 
galaxies. Dust evolution models seek to explain and 
quantitatively reproduce these variations. Most of these 
models treat the dust evolution problem averaged over a 
galaxy, and follow the basic template set by Audouze 
& Tinsley (1976) in which different species (e.g. hy¬ 
drogen or a specific metal) travels between well-mixed 
states (e.g. gas-phase outside of a galaxy or solid-phase 
in the dense molecular medium in a galaxy) according 
to a set of differential equations. The current genera¬ 
tion of dust evolution models (e.g., Dwek & Cherchneff 
2011; Zhukovska 2014; Feldmann 2014) include a more 
comprehensive and relatively well-constrained set of dust 
production and destruction mechanisms than the origi¬ 
nal, but are still focused on reproducing galaxy-averaged 
properties. As the number of destruction and production 
parameters that need tuning is quite large, there is more 
than one way to correctly reproduce a galaxy-averaged 
dust-to-gas ratio. The current-generation models listed 
above can all reproduce trends that have been observed 
between the metallicity and GDR of a galaxy while dis¬ 
agreeing on what the most important mechanism for dust 
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evidence of dust destruction and/or growth in the 
Magellanic Clouds, ISM: abundances, ISM: dust 

evolution is. One possible way to break this degeneracy 
is to increase the number of constraining observations by 
attempting to reproduce GDRs averaged over different 
phases of the interstellar medium (ISM) of each galaxy. 

Computing an ISM-phase-averaged GDRs over a 
galaxy requires taking a comprehensive inventory of gas 
and dust over multiple phases, which requires observa¬ 
tions of different gas mass tracers with sufficiently high 
spatial resolution to resolve the dense molecular phase. 
Some of the most complete necessary datasets are avail¬ 
able for the Large and Small Magellanic Clouds (LMC 
and SMC, together - the MCs), which are two nearby, 
sub-solar metallicity (1/2 and 1/5 of the solar metallic¬ 
ity, respectively) dwarf galaxies. The dust content of 
the MCs has been extensively studied as part of the 
SAGE (Meixner et al. 2006; Gordon et al. 2011) and 
HERITAGE (Meixner et al. 2013) programs. As part 
of these programs, dust mass (Gordon et al. 2014a) and 
GDR (Roman-Duval et al. 2014) maps of the LMC and 
SMC have been made. 

When these GDR maps are binned by ISM phase, one 
finds different GDRs in diffuse neutral gas and dense 
molecular gas. However, the difference between these 
GDR values is not quantitatively accurate. Due to degen¬ 
eracies that are explored in detail in Roman-Duval et al. 
(2014), it is not currently possible to accurately compute 
dust masses in dense molecular gas from dust emission. 
One way to resolve some of these degeneracies is to mea¬ 
sure GDRs in various environments in the LMC and SMC 
using a method whose systematic uncertainties are differ¬ 
ent from those of the dust emission method. One of the 
aims of this study is to measure GDRs in the LMC and 
SMC using elemental depletions derived from ultraviolet 
(UV) absorption spectroscopy. A depletion is the differ¬ 
ence between an element’s gas-phase abundance and its 
intrinsic (i.e. combined gas-phase and solid-phase) abun- 
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dance. If one makes the assumption that the missing 
amount is entirely in dust, which is a reasonable assump¬ 
tion in the molecule-poor diffuse neutral medium (DNM) 
through which UV spectroscopy is possible, depletions 
can be converted to solid-phase abundances. Combining 
the solid-phase abundances of the main constituents of 
dust yields a GDR. 

Almost all depletion studies to date have focused on 
the DNM of the MW. The two main astrophysical re¬ 
sults of the ensemble of MW depletion studies are that 
in a single location, elements with higher condensation 
temperatures tend to be more depleted (Field 1974) and 
that the depletions of every element track each other and 
the volume density of the gas they are associated with in 
a consistent and continuous way (Jenkins 2009). The lat¬ 
ter result also implies that the gas volume density and 
GDR in a small amount of DNM are correlated. The 
SMC is the next-best studied galaxy after the MW, with 
four sightlines along which some or all of Mg, Si, and Fe 
have been measured (Welty et al. 1997, 2001; Sofia et al. 
2006). These four measurements hint at possible differ¬ 
ences between the composition of dust in the SMC and 
MW, but are inconclusive. 

No other system has depletion measurements of mul¬ 
tiple important dust constituents along more than two 
sightlines. There have been several one-to-two sightline 
studies of local galaxies (James et al. 2014 and references 
therein) and many single sightline studies of damped 
Lyman-a systems (DLAs) (Rafelski et al. 2012 and refer¬ 
ences therein). Because intrinsic elemental abundances 
in most local galaxies and all DLAs are not available, 
these studies have relied on MW depletion patterns to 
interpret their observed gas-phase abundances. The 
present study will provide depletion patterns for two 
more galaxies at sub-solar metallicities, which may be 
more appropriate contexts for the interpretation of gas- 
phase abundances in low-metallicity systems. 

In this paper, we present new measurements of silicon 
(Si II), phosphorus (P II), chromium (Cr II), iron (Fe 
II), and zinc (Zn II) column densities in the LMC and 
SMC. These measurements are made using 16 new and 
84 archival spectra. Our sample, observations, and data 
reduction are described in Section 2. Our data analysis, 
which involves measuring column densities and comput¬ 
ing gas-phase abundances, depletions, solid-phase abun¬ 
dances, and GDRs, is described in Sections 3 and 4. A 
discussion of our results and a comparison of gas-phase 
abundances in the MCs and DLAs is presented in Section 
5. Our results are briefly summarized in Section 6. 

2. OBSERVATIONS AND ARCHIVAL DATA 

The spectra presented in this study are a combina¬ 
tion of new observations from the Cosmic Origins Spec¬ 
trograph (COS; Green et al. 2012) on the Hubble Space 
Telescope (HST) along 16 sight lines and supplementary 
archival spectra for these 16 and 69 additional sight lines 
from the Far Ultraviolet Spectroscopic Explorer {FUSE 
Sahnow et al. 2000). Figure 1 shows the positions of 
these sight lines in the galaxies. 

2.1. COS Observations 

The COS spectra of the Magellanic Cloud sight lines 
presented here were obtained between the dates of 2012 
October 12 and 2013 August 19 as part of HST program 


number 13004. All successfully acquired targets were 
faint enough to use the primary science aperture, and 
were dithered with all 4 FP-POS positions to improve the 
limiting signal-to-noise ratio (SNR). During the observa¬ 
tion of target AzV 388, HST was unable to acquire the 
guide star, resulting in a mis-pointing and a failed obser¬ 
vation. Observations were conducted with the near-UV 
(NUV) G185M grating at central wavelengths of 1921 
and 1953 A and a spectral resolution of 17 km sec“^. 

2.2. COS Targets 

Table 2 lists the 8 SMC and 8 LMC stars along the lines 
of sight targeted in our new COS observations with their 
sky coordinates, spectral types, V, E(B-V), and HI and 
H 2 column densities as reported in Welty et al. (2012). 
These targets will be referred to as the COS sample. 

The COS sight lines were chosen from a set of 285 Mag¬ 
ellanic Cloud targets with HI and/or H 2 column den¬ 
sities measured from archival HST and FUSE spectra 
presented in Welty et al. (2012). In order to explore the 
dependence of depletion on the column density {N) and 
phase of gas, sight lines were chosen to most completely 
cover the 27V(H2)/(iV(HI)-L2A(H2)) vs. A(HI)-L 2 Af(H 2 ) 
parameter space. The targets sample a range of total 
neutral hydrogen column densities from iV(Hl) -|-2N(H2) 
of ~ 10^° ® to ^ 10^^ cm“^ and molecular hydrogen frac¬ 
tion, 2iV(H2)/(A^(HI)-h2iV(H2)), of ~ 50% to < 1%. 
Few sight lines with high molecular fractions exist in the 
original sample of 285, so our COS sample is dominated 
(14/16) by those with molecular fractions of < 10%. 

2.3. FUSE Targets 

Data for the 69 supplementary sightlines, of which 32 
are towards the LMC and 37 are towards the SMC, were 
downloaded from the FUSE MC Legacy Project archive 
(Blair et al. 2009) and analyzed with no further pro¬ 
cessing. The FUSE MC Legacy Project archive contains 
spectra towards a total of 287 stars. From these, we se¬ 
lected stars whose continuum SNRs in the region from 
1140 to 1155 A were greater than 5 and towards which 
Welty et al. (2012) had detected H I and H 2 . These 
targets will be referred to as the FUSE-only sample. 

2.3.1. Ancillary literature column densities 

Two of the LMC and four of the SMC targets in the 
FUSE-only sample have been observed and analyzed be¬ 
fore (Roth & Blades 1997; Welty et al. 2001; Sofia et al. 
2006) . All of these targets have Zn and Cr measurements 
and three of the four SMC targets have Si measurements. 
Spectra towards Sk 155 and Sk 108 has been separately 
analyzed by both Welty et al. (2001) and Sofia et al. 
(2006). We adopt the values of the latter when they are 
given and the former otherwise. We note that this work 
is the first to use these measurements to compute deple¬ 
tions relative to hydrogen; the original authors did not 
have access to the hydrogen absorption measurements of 
Welty et al. (2012). 

2.4. Data Reduction 

Figures 2 and 3 show spectra towards SK —65 22 in 
the LMC and SK 116 in the SMC. Labels and vertical 
lines indicate the typical wavelengths of absorption lines 
of interest; FUSE and COS line spread functions (LSFs) 
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TABLE 1 

Locations, stellar parameters, and ISM parameters of the COS sample 


Target 

Galaxy 

RA (J2000) 

Dec (J2000) 

Spectral type 

V 

E(B-V) 

log(iV(H I)) 

log(A(H 2 )) 

SK 9 

SMC 

11.6360 

-73.1016 

05 V 

13.55 

0.15 

21.76 

17.03 

AzV 47 

SMC 

12.2145 

-73.4329 

08 III 

13.44 

0.13 

21.32 

18.54 

AzV 95 

SMC 

12.8400 

-72.7375 

07 III 

13.78 

0.14 

21.49 

19.40 

AzV 238 

SMC 

14.9813 

-72.2272 

09.5 III 

13.64 

0.13 

21.41 

15.95 

AzV 327 

SMC 

15.7939 

-72.0373 

09.5 Il-Ibw 

13.03 

0.13 

20.93 

14.79 

SK 116 

SMC 

16.2323 

-72.7800 

09 labw 

12.59 

0.13 

21.57 

18.53 

SK 143 

SMC 

17.7324 

-72.7156 

09.7 Ib 

12.83 

0.36 

21.00 

20.93 

AzV 476 

SMC 

18.4269 

-73.2915 

06.5 V 

13.52 

0.23 

21.85 

20.95 

SK -65 22 

LMC 

75.3462 

-65.8759 

06 Iaf+ 

12.08 

0.11 

20.58 

14.93 

SK -67 5 

LMC 

72.5789 

-67.6606 

09.7 Ib 

11.34 

0.14 

21.00 

19.46 

SK -71 45 

LMC 

82.8154 

-71.0695 

04-5 Ill(f) 

11.54 

0.16 

21.09 

18.63 

SK -66 172 

LMC 

84.2725 

-66.3598 

02 III 

13.13 

0.18 

21.25 

18.21 

SK -70 115 

LMC 

87.2069 

-70.0661 

06.5 laf 

12.24 

0.20 

21.30 

19.94 

BI 173 

LMC 

81.7914 

-69.1324 

08.5 11(f) 

12.96 

0.15 

21.34 

15.64 

BI 42 

LMC 

74.2536 

-66.4070 

08 V 

12.95 

0.19 

21.51 

18.69 

SK -68 135 

LMC 

84.4547 

-68.9171 

ON9.7 la-l- 

11.36 

0.26 

21.60 

19.87 


Note. — All values from Welty et al. (2012) and references therein. 
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Fig. 1.— HI column density maps of the LMC (left) and SMC (right). The positions of the 42 LMC and 46 SMC archival FUSE sightlines 
used in this study are marked with dots. The 8 LMC and 8 SMC sightlines that also have new COS observations are emphasized with 
large, open circles. 


are shown at the bottom left of each panel (Kruk et al. 
2002; Ghavamian et al. 2009). The relative velocities 
assumed for LMC and SMC absorption are based on the 
ranges given in Briins et al. (2005). Most of the other 
spectra in this study resemble the ones described above. 

The continua of SK —68 135 and SK —65 22 con¬ 
tain broad absorption and emission lines, including some 
which overlap with absorption lines of interest. In some 
of the FUSE spectra, the absolute flux level between 
1140 and 1150 A is decreased relative to the flux level at 
nearby wavelengths due to a shadow from a repeller grid 
falling on part of the FUSE detector. All of our Ee II lines 
fall within this region; spectra where the shadow is espe¬ 
cially strong are only marginally usable and have been 
excluded. There are absorption lines due to gas between 
the MW and LMC in half of the observations towards the 
LMC. This gas has relative velocities which range from 


100 to 180 km sec~^, has been observed before (Lehner 
et al. 2009), and can usually be distinguished from MW 
and LMC absorption. 

The resolutions of EUSE and COS are ~ 13 km sec“^ 
(Kruk et al. 2002) and ~ 15 km sec~^ (Ghavamian et al. 
2009), respectively. Isolated lines in higher resolution 
(^ 3 km sec” ) spectra of the diffuse neutral medium 
in the MCs have widths of order a km sec” (Roth & 
Blades 1997; Welty et al. 1997), implying that all of the 
absorption we see is almost certainly unresolved. 

3. QUANTITATIVE ANALYSIS: GAS-PHASE 

We compute four measurements of gas-phase metal 
content, which are listed below in order of increasing 
distance from the data. From the spectra, we measure 
column densities for P II, Zn II, Si II, Gr II and Fe II. Nor¬ 
malizing these column densities by the hydrogen column 
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Fig. 2. — FUSE (top two panels) and COS (bottom three panels) observations of the LMC 06 Iaf+ star Sk -65 22. Labels show the 
locations of ISM absorption features of interest. The instrumental LSF for each wavelength range is shown at the bottom left of each panel. 
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Fig. 3.— FUSE (top two panels) and COS (bottom three panels) observations of the SMC 09 labw star Sk 116. Labels show the locations 
of ISM absorption features of interest. The instrumental LSF for each wavelength range is shown at the bottom left of each panel. 
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density gives us ion abundances. If we had measurements 
of a single element in more than one ionization state, we 
would compute an ionization correction for each sight¬ 
line. Instead, we check that the ionization corrections 
for most of our sightlines should be negligible and adopt 
the ion abundances and elemental abundances. Normal¬ 
izing these elemental abundances by reference, or total 
gas- and solid-phase, abundances gives us elemental de¬ 
pletions. Finally, we analyze the distribution of deple¬ 
tions in each galaxy using the formalism from Jenkins 
(2009). This formalism gives us a concise summary of 
an element’s depletion trends each galaxy and allows us 
to impute partially missing data, such as Zn, Si, and Cr 
depletions for sightlines in the FUSE-only sample. Be¬ 
low, we describe each of these steps in moderate detail; 
the detailed mechanics of steps with relatively complex 
implementations are described in the appendices. 

3 . 1 . Ion column densities 

The analysis of an absorption spectrum can generally 
be divided into two steps - continuum fitting and column 
density recovery. In a typical continuum fitting proce¬ 
dure, one interpolates over ISM absorption using neigh¬ 
boring parts of the spectrum as a guide. We do this inter¬ 
polation using Gaussian process regression (Rasmussen 
& Williams 2005). A Gaussian process is a probability 
distribution over functions. Regression is a procedure for 
choosing a function that has a high probability of having 
produced a set of observations. Gaussian process regres¬ 
sion, then, is a procedure for choosing a function that, 
on the one hand, has a high probability of having pro¬ 
duced the observed points and, on the other, has a high 
probability according to the specified Gaussian process. 
The procedure also produces an estimate of the func¬ 
tion’s pointwise uncertainty. This estimate includes cor¬ 
relations between neighboring values. Our chosen Gaus¬ 
sian process favors functions that resemble the sum of 
a second-order polynomial and a smooth (i.e. infinitely 
differentiable), zero-mean perturbation. For an example 
of what these continua look like over short wavelength 
ranges, see the dashed black line in Figure 4. For more 
details on this procedure, see Appendix A. 

The next step, column density recovery, requires a pro¬ 
cedure for converting the observed amount of absorption 
into a column density or, if different parts of the absorp¬ 
tion are caused by different absorbers, column densities. 
In our spectra, the absorption line corresponding to each 
transition appears at several Doppler velocities, each cor¬ 
responding to a different line-of-sight object. In Figure 
4, each absorption line appears near 0 km sec“^, which 
corresponds to gas in the MW, and 130 km sec“^, which 
corresponds to gas in the SMG (Briins et al. 2005). To¬ 
wards stars in the LMG, we see absorption at 0, 150, and 
250 km sec~^ due to gas in the MW, intervening high 
velocity clouds (Lehner et al. 2009), and and the LMC. 
These shifts cause velocity components of some absorp¬ 
tion lines, in particular the Gr II and Zn II lines near 
2062 A and the three Fe II lines near 1143 A, to overlap. 
This can be seen in the second and third panels from the 
bottom of the right column of Figure 4 in which each 
absorption line has two broad velocity components, but 
appears to have four. 

To resolve the resulting confusion, we turn to a vari¬ 


ation on Voigt profile fitting. In standard Voigt profile 
fitting, one models the observed absorption as coming 
from a series of coherent clouds. Each cloud has a cen¬ 
tral velocity and (Gaussian) velocity width and contains 
some amount of each species under investigation. One 
optimizes the central velocity, width, and amount of each 
species in each cloud by generating absorption spectra, 
comparing them to the observations, and minimizing the 
difference. The number of clouds is chosen in order to get 
an acceptable fit according to some criterion. We avoid 
choosing a set number of clouds and instead numerically 
integrate out the dependence on the number of coherent 
clouds using a transdimensional (i.e. operating over a 
parameter space of varying size) Markov Ghain Monte 
Garlo (MGMG) procedure (Green 1995). This procedure 
allows us to include the effects of different possible cloud 
decompositions in our quoted uncertainties. Additional 
information about our Voigt profile fitting implementa¬ 
tion can be found in Appendix A. 

We have used these two procedures to measure ion col¬ 
umn densities. These are plotted in Figure 5 and listed in 
Table 3 (COS sample) and Table 3 (FUSE-only sample). 

3 . 2 . Elemental abundances and ionization corrections 

Next, we use the measured ion column densities to 
compute gas-phase elemental abundances. From our 
ion measurements and literature atomic and molecular 
hydrogen column densities (Welty et al. 2012), we can 
directly compute gas-phase ion abundances. Along a 
sightline with Zn II column density iV(Znll), atomic hy¬ 
drogen column density V(HI), and molecular hydrogen 
column density iV(H 2 ), the gas-phase Zn II abundance 
egas(Znll) is 

/ V(Znll) \ 

Sgasi n )= Ogio 2 iV(H 2 )) 

Because this approximation does not account for every 
possible gas-phase form of zinc and hydrogen, the Zn II 
abundance may be different from the true gas-phase Zn 
abundance £gas(Zn); the difference between the ion and 
element abundances is called the ionization correction: 

IG(Znll) = SgaA^n) - Ega^Znll). (2) 

The magnitude of the ionization correction depends on 
the element’s ionization potentials and the ISM phases 
present along the line of sight (Sembach et al. 2000). 
Based on models of the ionization structure of clouds 
with H I and H 2 column densities similar to those mea¬ 
sured along our sightlines, we expect the ionization cor¬ 
rections to be negligible. 

The ionization correction for the Zn II abundance is 
small when, on the one hand, most of the Zn that is 
spatially coincident with H I and H 2 is Zn II and, on the 
other hand, most of the Zn that is not spatially coincident 
with H I and H 2 is not Zn II. If the main ionizing source is 
the background star towards which we are observing and 
its neighbors, then both conditions for a small ionization 
correction will apply. Both conditions require radiation 
fields that should affect a qualitatively small fraction of 
the total gas. Because the ionization potential from H I 
to H II is lower than the ionization potential from Zn II 
to Zn III, the only way to miss some of the spatially co¬ 
incident Zn is to attenuate the radiation field above the 
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Fig. 4.— Profile fit of the observation towards Sk 116 in the SMC for the different ions. Listed next to the ions is the rest frame {v = 0) 
wavelength of the absorption line. Circles represent the data; light and dark grey regions are the 68 and 95% credible regions about the 
median fit; the dashed line is the median continuum level 
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TABLE 2 

Logarithmic ion column densities towards LMC and SMC targets in the COS sample 


Target 

Galaxy 

logio(iV(PII)) 

[cm 

logio(iV(ZnII)) 

[cm 

logio(Af(SiII)) 

[cm 

logio(iV(CrII)) 

[cm 

logio(iV(FeII)) 

[cm 

logic (-^h) 
[cm“^| 

AzV 327 

SMC 

13.30 ± 0.13 

12.75 ±0.19 

15.69 ±0.04 

11.74 ± 1.02 

14.67 ± 0.02 

20.93 

SK 143 

SMC 

12.33 ±0.99 

12.57 ± 0.33 

15.10 ±0.11 

11.49 ±0.97 

14.65 ±0.03 

21.43 

AzV 47 

SMC 

13.77 ± 0.03 

13.08 ±0.03 

15.74 ±0.02 

13.66 ±0.03 

15.43 ±0.02 

21.32 

AzV 238 

SMC 

13.65 ±0.04 

13.05 ±0.04 

15.73 ±0.02 

13.31 ± 0.07 

14.91 ± 0.02 

21.41 

AzV 95 

SMC 

13.77 ± 0.01 

13.15 ±0.03 

15.86 ±0.02 

13.61 ± 0.04 

15.30 ±0.01 

21.50 

SK 116 

SMC 

13.83 ±0.04 

13.15 ±0.02 

15.95 ±0.02 

13.66 ±0.02 

15.36 ±0.03 

21.57 

SK 9 

SMC 

13.09 ± 1.38 

13.33 ±0.02 

16.04 ±0.03 

13.70 ± 0.04 

13.91 ± 2.10 

21.76 

AzV 476 

SMC 

14.02 ±0.04 

13.38 ±0.04 

16.04 ± 0.04 

13.75 ±0.04 

15.41 ± 0.04 

21.95 

SK-65 22 

LMC 

13.67 ± 0.01 

12.69 ±0.12 

15.57 ± 0.02 

13.10 ± 0.04 

14.83 ±0.01 

20.58 

SK-67 5 

LMC 

13.73 ±0.02 

13.03 ±0.02 

15.48 ±0.02 

13.22 ±0.03 

14.56 ±0.02 

21.02 

SK-71 45 

LMC 

14.02 ±0.01 

13.17 ± 0.04 

16.01 ± 0.01 

13.61 ± 0.03 

15.25 ±0.01 

21.09 

SK-66 172 

LMC 

13.95 ±0.03 

13.07 ± 0.04 

15.65 ±0.02 

13.10± 0.10 

14.68 ±0.04 

21.25 

SK-70 115 

LMC 

13.97 ± 0.02 

13.16 ±0.03 

15.90 ±0.01 

13.67 ± 0.02 

15.39 ±0.02 

21.34 

BI 173 

LMC 

13.70 ± 0.02 

13.13 ±0.02 

15.73 ±0.02 

13.54 ± 0.03 

15.04 ± 0.02 

21.34 

BI 42 

LMC 

13.90 ± 0.11 

13.24 ± 0.02 

15.91 ± 0.09 

13.52 ±0.04 

15.35 ±0.08 

21.51 

SK-68 135 

LMC 

14.06 ±0.05 

13.36 ±0.06 

15.90 ±0.06 

13.76 ±0.02 

15.45 ±0.03 

21.62 
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Fig. 5. — Metal ion column density plotted against the total 
column density of atomic and molecular hydrogen. LMC sightlines 
are denoted in red, SMC sightlines are in blue. Sightlines without 
COS data are shown without errorbars. 


TABLE 3 

Logarithmic ion column densities towards LMC and SMC 

TARGETS IN THE FUSE SAMPLE 


Target 

Galaxy 

logio(Af(PII)) 

[cm 

logio(Af(FeII)) 

[cm 

logio(^H) 

[cm-2[ 

AzV 6 

SMC 

14.09 ±0.12 

15.40 ± 0.05 

21.54 

SK 10 

SMC 

13.95 ± 0.06 

15.50 ±0.02 

21.58 

SK 15 

SMC 

14.25 ± 0.35 

15.61 ± 0.13 

21.80 

SK 18 

SMC 

14.09 ± 0.06 

15.51 ± 0.04 

21.77 

SK 34 

SMC 

14.29 ±0.42 

15.39 ±0.05 

21.59 


Note. — Hydrogen column densities are from Welty et al. (2012) 
and references therein. 


Table 3 is published in its entirety in the electronic edition of the 
Astrophysical Journal. A portion is shown here for guidance regarding 
its form and content. 


ionization potential from Zn 1 to Zn II. Since we can see 
non-negligible fluxes at and below the wavelength cor¬ 
responding to this ionization potential, 1320 A, the gas 
along the line of sight cannot contain significant amounts 
of Zn I. To violate the second condition, that Zn II be 
present without H I or H 2 , we need a radiation field that 
is capable of ionizing H I but not Zn II. This requires 
photons with wavelengths between 690 A and 912 A, 
which corresponds to a fairly narrow shell within an H II 
region. The ionization potentials to and from the other 
ions in this study imply similar qualitative limits on the 
corresponding ionization corrections. 

We can make these qualitative limits quantitative us¬ 


ing models of clouds illuminated by UV-bright stars. We 
do not know the volume density of the gas along each 
sightline, so we have cloud models with constant vol¬ 
ume densities of 0.1, 1, 10, and 100 particles per cm^. 
The spectral types of our background stars range from 
BO to 02, so we run each cloud model with a BO star 
and with an 02 star in order to bracket the range of 
possible ionization corrections; the stellar spectra come 
from Lanz & Hubeny (2003). Stellar spectra, gas cool¬ 
ing, and gas shielding all depend on the metallicity, so we 
run each star-cloud pair at the LMC and SMC metallic- 
ities. These and other model parameters are listed and 
explained in Appendix B. For each combination of spec¬ 
tral type, volume density, and metallicity, we use Cloudy 
(version 13.02, Ferland et al. 2013) to calculate the cloud 
ionization structure. 

From the output of these calculations, we derive ion¬ 
ization corrections and molecular hydrogen fractions as a 
function of H I and H 2 column density, where the molec¬ 
ular hydrogen fraction /h 2 is 


2iV(H2) 

2Af (H 2 )-biV(HI)’ 


(3) 


Ionization corrections for P II and Zn II are shown in Fig¬ 
ure 6, along with the positions of our observations in the 
hydrogen column-molecular fraction plane. These cor¬ 
rections are for an 02 star in the LMC. Almost all of the 
stars towards which we observe have later spectral types 
and, correspondingly, lower-magnitude P II corrections. 
The P II ionization corrections for this specihc combi¬ 
nation of metallicity and illuminating star have larger 
absolute values than the ionization corrections for any 
other ion in any other model. In all cases other than P 
II at the LMC metallicity illuminated by an 02 star, the 
ionization corrections are similar in sign and magnitude 
to those of Zn II in Figure 6, whose largest magnitude in 
the part of the diagram occupied by our observations is 
less than 0.1 decades. To within our typical uncertain¬ 
ties, which are about 0.1 decades, our ion abundances 
should be equal to the gas-phase elemental abundances. 


3.3. Depletions 

The depletion (5(X) of an element X is the fraction of 
the element’s assumed reference, or gas- and solid-phase 
total, abundance e^ef that we observe in the gas phase: 


(5(X) = logio 



— e(X)gas e(X)ref, 


(4) 

(5) 


where Cgas is the observed gas-phase abundance. The 
more negative 5 is, the more of the element is missing, 
presumably in the solid phase. This makes 5 a proxy for 
the gas-to-dust ratio. 

In order to compute depletions, we need reference 
abundances. The composition of a star’s photosphere 
is a superposition of the ISM composition at the star’s 
formation time with the effects of various enrichment 
processes. This makes recently-formed stars that have 
not yet undergone self-enrichment good proxies for the 
present-day ISM composition. 
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Fig. 6.— Theoretical ionization corrections (colors) for P II (left) and Zn II (right) as a function of the non-ionized hydrogen column 
density and the molecular hydrogen fraction. Black points mark the hydrogen column densities and molecular hydrogen fractions along 
sightlines with FUSE (left panel) and COS (right panel) data. Each track corresponds to a single total hydrogen volume density, which is 
marked at the track’s base in left panel. These corrections would apply to a sightline towards an 02-type star in the LMC. The corrections 
for every element in our study other than P II are similar to those of Zn II. See Subsection 3.2 for a description of our photoionization 
calculations. 


TABLE 4 

PhOTOSPHERIC abundances in the LMC AND SMC 


Element 

LMC abundance 
[cm“^] 

SMC abundance 
[cm“^] 

References 

c 

7.94 ±0.10 ±0.05 

7.52 ±0.10 ±0.04 

1, 2, 3, 5, 6, 

7, 8, 10, 11 

o 

8.50 ±0.11 ±0.03 

8.14 ±0.08 ±0.04 

1, 2, 3, 5, 6, 

9, 10, 12 

Mg 

7.26 ±0.08 ±0.04 

6.88 ± 0.06 ± 0.03 

1, 2, 4, 5, 6, 

9, 10, 11, 12 

Si 

7.35 ±0.10 ±0.02 

6.96 ±0.07 ±0.09 

1, 2, 4, 5, 6, 

9, 10, 12 

P 

5.1 

4.7 

0 

Cr 

5.37 ± 0.07 ± 0.03 

4.92 ±0.10 ±0.05 

1, 2, 4, 9, 11, 12 

Fe 

7.32 ± 0.08 ± 0.03 

6.89 ± 0.08 ± 0.03 

1, 2, 4, 6, 9, 11, 12 

Ni 

5.92 ±0.07 ±0.03 

5.52 ±0.18 ±0.07 

1, 2, 4, 9, 11 

Zn 

4.31 ± 0.15 ± 0.15 

4.02 ±0.20 ±0.14 

8, 9, 11 


Note. — Abundances are given in the form log 2 ^Q(X/H) + 12. The first 
uncertainty is the standard deviation of the abundance value; the second is the 
standard deviation of the star-to-star dispersion. 

References. (0) Solar (Grevesse et al. 2010), (1) Chekhonadskikh 

(2012), (2) Hill et al. (1995) (3) Hill et al. (1997) (4) Hill (1997) (5) Hunter 
et al. (2007) (6) Korn et al. (2002) (7) Korn et al. (2005) (8) Luck & Lambert 
(1992) (9) Luck et al. (1998), (10) Rolleston et al. (2003), (11) Russell & Bessell 
(1989) (12) Venn (1999) 


There have been a number of spectroscopic measure¬ 
ments of the composition of luminous stars in both MCs. 
Because no individual study includes all of the elements 
we are interested in, we have to combine their mea¬ 
surements. We believe that we cannot simply average 
measurements across the studies because of significant 
between-study differences, repeat observations and anal¬ 
yses of the same stars by multiple studies, and miss¬ 
ing uncertainty information. Instead, we partially pool 
measurements across studies using a multilevel linear 
model (Browne & Draper 2006). The model is described 
in detail in Appendix C. Our adopted reference abun¬ 
dances, abundance uncertainties, and estimated intrinsic 
(i.e. location-to-location) abundance variances are listed 
in Table 3.3. Figure 7 shows the depletions we obtain by 
assuming these reference abundances. 


In the next section, we analyze how depletions vary 
from sightline to sightline. 

3.4. Linear depletion relations 

We would like to use the depletions derived in the pre¬ 
vious section to estimate ranges of solid-phase elemental 
abundances in the MCs and compare them with abun¬ 
dances in the MW. For these estimates to be robust, we 
need a representative sample of each MC’s diffuse neutral 
medium. The difference between the lowest and highest 
depletions of Fe and P, which we measure using FUSE 
spectra, is approximately equal to the same difference 
in the MW, suggesting that the FUSE sample is repre¬ 
sentative of the same range of diffuse neutral medium 
conditions as is covered by MW datasets. 

In order to make the same depletion range estimates 
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Fig. 7.— Elemental depletions relative to the ISM abundances 
of the LMC and SMC plotted against the total column density of 
atomic and molecular hydrogen. Depletions towards the LMC and 
SMC are shown in red and blue. 


TABLE 5 

Abundance uncertainties and depletion parameters 



LMC 

SMC 

P: 



5(P)o 

0.30 ±0.03 

-0.30 ± 0.04 

Ap 

-0.71 ± 0.07 

-0.83 ±0.07 

Zn: (Je 

0.15 

0.2 

<5(Zn)o 

-0.08 ±0.06 

-0.04 ± 0.06 

^Zn 

-0.60 ±0.06 

-0.61 ± 0.06 

Si: (Tg 

0.1 

0.07 

<5(Si)o 

-0.32 ±0.03 

-0.01 ± 0.08 

Asi 

-0.92 ±0.08 

-1.03 ±0.15 

Cr: 

0.07 

0.1 

<5(Cr)o 

-0.77 ± 0.04 

-0.36 ±0.06 

Acr 

-0.71 ± 0.09 

-0.97 ± 0.11 

Fe: (Tg 

0.08 

0.08 

5(Fe)o 

-0.96 ±0.02 

-0.57 ± 0.05 

Ape 

-1.00 ±0.06 

-1.13 ±0.07 


Notes. We describe the abundance uncertainties Ce in Section 3.3 
and the depletion parameters (5(X)o and Ax in Section 3.4. 

for the elements that are only available in the (unrepre- 
sentative) COS sample, we need a way of combining its 
elemental information with the FUSE sample’s popula¬ 
tion information. For this purpose, we use the Jenk¬ 
ins depletion parametrization (Jenkins 2009). In this 
parametrization, the observed depletion (5(A)i of element 
X along a sightline i is given by 

( 6 ) 


where i5(X)o and Ax are the minimum and range of de¬ 
pletions of element X in the diffuse neutral medium, and 
F* is a parameter representing the overall depletion level 
of sightline i. F* is kept fixed for all elements along sight¬ 
line i and J(X)o and Ax are kept fixed for an element X 
across all sightlines. 

This model is underdetermined; we can rescale or shift 
combinations of parameters without changing the obser- 
ables S{X)i. We can scale a set of slope parameters Ax 
by a constant C and scale the F* values by 1/C. At fixed 
Ax, we can add a constant C to i5(X)o and subtract a fac¬ 
tor of C/Ax from the F*. We fix these two degeneracies 
and connect our F* values with those of Jenkins (2009) 
by imposing the MW values as priors on i5(Zn) and Azn- 
We chose this prior because its implied F* zero point lo¬ 
cation has a clear interpretation - F* is zero when the 
depletion of Zn is (up to uncertainty) zero. Since Zn is a 
volatile element in the MW, meaning that its minimum 
depletion in the DNM is zero, we expect it to continue 
to be volatile in more dust- and metal-poor galaxies such 
as the MCs. 

Our priors on the remaining parameters are uninfor¬ 
mative. We impose broad uniform priors on the other 
depletion zero points and slopes, a uniform prior with 
bounds and E^ax s-^id uniform priors 

with bounds (—3, 0.5) and (0.5, 3) on, respectively, FF^ 
and The F* bounds are included as parameters 

because we do not a priori know how much of the pos¬ 
sible F* range our observations will fill. For both MCs, 
the F* ranges ended up being approximately (0,1). 

We implemented and generated samples from this 
model using the Bayesian statistical analysis module 
PyMC (Path et al. 2010). Our fits to the LMC and SMC 
depletions are shown in Figures 8 and 9. For each el¬ 
ement, the reference abundance uncertainties, depletion 
zero points (5(X)o, and depletion slopes Ax are listed in 
Table 5. The depletion zero points of the refractory ele¬ 
ments Si, Cr, and Fe are significantly less negative in the 
SMC than they are in MW or LMC. The minimnm Si de¬ 
pletion is particularly notable, since it is consistent with 
0. However, if we include the uncertainty of the reference 
Si abundance, the minimum Si depletion is also within 
1(7 of 0.13 in logarithmic units, which corresponds to 1/4 
of the Si being out of the gas phase. The depletion slope 
of Cr in the LMC and SMC is similar and signficantly 
smaller than in the MW. The other depletion slopes do 
not change significantly between the three galaxies. In 
the common interpretation of depletion variations within 
a galaxy, depletion zero points are associated with the 
composition of the refractory dust cores while depletion 
slopes are associated with the rate at which dust mantles 
either grow or are destroyed. By this interpretation, the 
fractional gas-to-solid transfer rate is roughly constant 
between galaxies while the refractory core composition 
changes. We will discuss the implications of these zero 
point and slope changes further in Sections 5.3 and 5.4. 

4. QUANTITATIVE ANALYSIS: SOLID-PHASE 

In this section, we describe our derivation of solid- 
phase abundances and GDRs from measured gas-phase 
and assumed total ISM abundances. 

4.1. Solid-phase abundances 


6{X), = 5{X)o + AxF*, 
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Fig. 8.— Fits of the Jenkins parametrization to elemental depletions in the LMC ISM. Data points include our observations, our analysis 
of archival FUSE data, and measurements from Roth & Blades (1997). Shaded regions show 68% credible intervals about the median fit. 
The black line is the MW best fit from Jenkins (2009). 
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Fig. 9. — Fits of the Jenkins parametrization to elemental depletions in the SMC ISM. Data points include our observations, our analysis 
of archival FUSE data, and measurements from Roth & Blades (1997), Welty et al. (2001), and Sofia et al. (2006). Shaded regions show 
68% credible intervals about the median fit. The black line is the MW best fit from Jenkins (2009). 
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Fig. 10.— We show the surface mass GDR as a function of the 
over-all depletion level F *, which increases from low depletion in all 
elements {F* = 0) to high depletion in all elements {F* = 1). The 
MW, LMC, and SMC are denoted by black, red, and blue. The 
solid lines are derived by applying Equation 10 to each galaxy’s 
best-fit reference abundances and the MW’s best-fit depletion re¬ 
lations (Jenkins 2009). The shaded regions are derived by apply¬ 
ing Equation 10 to each galaxy’s reference abundances and deple¬ 
tion relations, both with uncertainties. The vertical extent of the 
shaded regions represent the 84% posterior credible interval of the 
GDR at each value of F*. 


We can reasonably assume that all of the metals that 
are missing from the gas-phase are in dust and can 
use this assumption to derive a solid-phase abundance 

^(^)dust- 

£(X)dust = logio (7) 

= logio(l0"^^^"^ , (8) 

where £(X)ref is the reference (i.e. theoretical gas- and 
solid-phase total) abundance, £(X)gas is the gas-phase 
abundance, and (5(X) is the depletion of the gas-phase 
abundance relative to the reference abundance. By com¬ 
bining the depletion-based formulation of Equation 8 
with the depletion parametrization described in Sec¬ 
tion 3.4, we can get a continuous empirical relation be¬ 
tween F* and the solid-phase abundance of an element. 
We explore solid-phase abundances and related quanti¬ 
ties further in Sections 4.2 and 4.3. 

4.2. GDRs 


We can combine all of the elements’ solid-phase abun¬ 
dances to compute a gas-to-dust mass ratio: 


1 

GDR 


-'dust 


^gas 




(9) 


where X ranges over all elements and mx is the atomic 
mass of element X. This expression can be rewritten in 
terms of depletions d(V): 

A_io^(x)y (10) 

^gas ^ 


In the MW, the only elements that contribute signifi¬ 
cantly to this sum are C, O, Mg, Si, Fe, and Ni. All other 
elements have low reference abundances and/or negligi¬ 
ble depletions. If we assume that the composition of dust 
in the MCs is even remotely like that of dust in the MW, 
we can restrict the sum in Equation 9 to those elements 
without missing much of the total dust mass. 

Since most of our sightlines do not have Si measure¬ 
ments and all of our sightlines do not have C, O, Mg, or 
Ni measurements, we compute GDRs using the linear de¬ 
pletion relations from Section 3.4. We use our LMC and 
SMC relations for Si and Fe and the MW relations from 
Jenkins (2009) for C, O, Mg and Ni. GDRs derived from 
these combination of depletion relations are shown as red 
(LMC) and blue (SMC) shaded regions in Figure 10. For 
comparison, we show the MW’s GDR in black and GDRs 
derived by assuming MW depletions at LMC and SMC 
reference abundances as solid red and blue lines. We will 
refer to these GDRs as abundance-scaled MW GDRs to 
differentiate them from the first set, which we will refer 
to as the LMC and SMC GDRs. 

We can make three main observations about this figure. 
First, regardless of which set of depletions we adopt, the 
GDR of all three galaxies changes by a factor of 2 from 
F* = 0 to F* = 1. This suggests that there is a signif¬ 
icant amount of dust destruction and/or growth in the 
DNM of the LMC, SMC, and MW. Second, the GDRs 
of the LMC and the high-depletion part of the SMC are 
approximately the same as the corresponding abundance- 
scaled MW GDRs. Third, the GDR of the low-depletion 
part of the SMC is higher than the SMC-abundance- 
scaled MW GDR. The second observation may be an 
artifact of our assumptions, as the MC depletion rela¬ 
tions that we do observe are always either as depleted or 
less depleted than the corresponding MW relations. This 
suggests that the depletions we are assuming for some or 
all of C, O, Mg, and Ni are too high and implies that 
even the GDRs which assume MC depletions for Si and 
Fe are most likely underestimates. The third observation 
suggests that the dust-to-metals ratio (DMR) decreases 
with decreasing metallicity. 

At and above the metallicity of the SMC, the (galaxy- 
averaged) DMR is constant while below the SMC metal¬ 
licity, the DMR ratio decreases with decreasing metal¬ 
licity (Remy-Ruyer et al. 2014). From the first observa¬ 
tion, we conclude that the DMR decreases by a factor of 
at least 2 within a galaxy. At medium to high over-all 
depletion (i.e. at F* = 1 or 0.5), the DMR does not 
change from the MW to the LMC to the SMC. However, 
at low over-all depletion (i.e. at F* =0), the DMR is 
constant from the MW to the LMC but decreases from 
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Fig. 11. — We show the solid-phase Fe abundance as a function 
of the solid-phase Si abundance. Data points from the SMC (blue) 
and LMC (red) are derived from reference abundances for the ap¬ 
propriate galaxy and gas-phase abundances using Equation 7, and 
data point error bars include uncertainties in both quantities. The 
MW region (black) is derived from solar abundances and MW de¬ 
pletions (Jenkins 2009) using Equation 8. The dashed grey lines 
mark levels of constant stoichiometric Fe to Si ratio. Within each 
galaxy, the over-all depletion level increases from left to right. 

the LMC to the SMC. This suggests that in diffuse, low- 
depletion gas, the DMR stops being constant at a higher 
metallicity than in more dense, higher-depletion gas. 

4.3. The solid-phase stoichiometry of silicon and iron 

As shown in the previous section, the GDR is not con¬ 
stant throughout the DNM of each galaxy. We now use 
depletion measurements to show that throughout each 
galaxy’s DNM, the dust stoichiometry, or its relative bulk 
composition, also varies. Using Equation 7, we have com¬ 
puted the solid-phase abundance of Si and Fe towards 
each sightline for which gas-phase column density mea¬ 
surements are available. These solid-phase abundances 
are shown in Figure 11, along with solid-phase abun¬ 
dances corresponding to the observed range of Si and 
Fe depletions in the MW (Jenkins 2009). The solid- 
phase Fe abundance is approximately constant across 
each galaxy, while the solid-phase Si abundance changes 
from the least- to the most-depleted sightlines by a factor 
of about 2. 

In all three galaxies, the Fe:Si abundance ratio in dust 
ranges from about 2:1 along the least-depleted sightlines 
to about 1:1 along the most-depleted sightlines. Since 
we are merely measuring the bulk dust composition, we 
cannot use these ratios to determine the composition 
of each dust subspecies. Nevertheless, the constant Fe 
abundance and variable Si abundance suggest that the 
process by which most of the Fe forms into dust is dif¬ 
ferent from the process by which at least half of the Si 
forms into dust. This suggests that the Fe-rich dust com¬ 
ponent is either particularly resilient, if most of the abun¬ 
dance variations are due to dust destruction, has a faster 



[X/H] 


Fig. 12. — Bars show the distribution of measured gas-phase P 
(dark color) and Zn (light color) abundances relative to each ele¬ 
ment’s solar abundance; vertical lines show the typically assumed 
metallicities of the LMC and SMC. P and Zn are relatively volatile 
elements. Their abundance at minimum depletion is a proxy for the 
ISM metallicity. Abundances greater than the mean MC metallici¬ 
ties (vertical lines) would indicate small-scale metallicity enhance¬ 
ments. We see no evidence for localized metal enrichment along 
our sightlines. 

growth rate, if most of the abundance variations are due 
to ISM dust growth, or some combination of the two. 
While this last fact has been known for MW dust for 
quite some time (see, e.g. Savage & Sembach 1996), we 
can now also confirm that it holds in other galaxies. 

5. DISCUSSION 

5.1. The ISM metallicity in the Magellanic Clouds 

Studies of field star abundances in the MCs (e.g. Car¬ 
rera et al. 2008; Cioni 2009) have found mild metallicity 
gradients as a function of galactocentric distance. The 
metallicities of stars of similar ages but different galacto¬ 
centric radii tend to be the same, suggesting that these 
gradients are due to star formation occuring at differ¬ 
ent times in different locations. If the observed gradi¬ 
ents were instead due to ISM inhomogeneities, stars that 
formed at the same time in different locations would have 
different metallicities. 

We can use the abundances of our more volatile ele¬ 
ments (P and Zn) to estimate the amount of localized 
metal enrichment towards some of our sightlines. At the 
lowest depletion levels, the gas-phase abundances of P 
and Zn should approach their total ISM abundances. An 
undepleted metal-enriched sightline would have a (local) 
metallicity greater than its corresponding galactic mean. 
Figure 12 shows the distribution of P and Zn abundances 
in the MCs relative to the solar abundance. While there 
are six apparently P-enriched and one apparently Zn- 
enriched sightlines across the pair of galaxies, this does 
not We do not find conclusive evidence of localized metal 
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Fig. 13. — Fe (from this work) and Ti (Welty & Crowther 2010) 
depletions in the LMC (red) and SMC (blue). For comparison, the 
linear depletion relation for Fe and Ti in the MW from Jenkins 
(2009) is shown in black. 

enrichment. 

Since even the relatively volatile elements in this study 
can be significantly depleted, we cannot exclude the pos¬ 
sibility that some of our sightlines are metal-poor. Our 
maximum measured gas-phase P and Zn abundances are 
consistent with the typically assumed MC metallicities 
of 1/2 and 1/5 times the solar metallicity for the LMC 
and SMC, respectively. 

5 . 2 . Previous depletion studies of the Magellanic Clouds 

Depletion studies of the MCs have been either medium¬ 
sized (< 20 depletion measurements) surveys of single el¬ 
ements or multi-element case studies of single sightlines. 
All of the targets in the multi-element case studies are 
part of our FUSE sample, and have been included in our 
analyses as described in Section 2.3. 

Welty & Crowther (2010) presented an optical survey 
of titanium (Ti) towards approximately twenty stars in 
either MC. Figure 13 shows a comparison of our Fe de¬ 
pletions and their Ti depletions for the sightlines which 
are included in both surveys. LMC and SMC sightlines 
are shown in red and blue. The black line represents the 
typical correlation of Fe and Ti depletions in the MW 
(Jenkins 2009). Across both galaxies, all but two of the 
targets are above this line, meaning that at fixed Fe de¬ 
pletion, MW sightlines tend to contain a smaller fraction 
of gas-phase Ti than MC sightlines. 

5 . 3 . The evolution of gas-phase abundances with 
metallicity 

Gas-phase abundance variations between pencil beam 
observations within a single galaxy can be succinctly 
summarized by the per-sightline depletion strength F* 



[Zn/H] 


Fig. 14.— Blue and red data points are abundance measurements 
for the SMC and LMC from our COS sample. The black line is 
based on a fit to gas-phase abundances in the DNM of the MW 
(Jenkins 2009). In the top panel, the red (LMC) and blue (SMC) 
lines are the MW line shifted to each galaxy’s metallicity. In the 
bottom panel, red and blue lines are the LMC and SMC’s linear 
depletion relations (see Sec. 3. 4). Since all of the LMC and SMC 
points in the top panel are above their corresponding lines, gas- 
phase abundances in the MCs differ from those of the MW by 
more than just the relative abundance differences. 

of the linear depletion relations of Section 3.4. Be¬ 
cause these within-galaxy variations are quite large, 
with depletion-driven gas-phase abundance differences 
between sightlines of up to of order a decade, galaxy-to- 
galaxy comparisons should be not of individual sightlines 
but of galaxy-wide depletion relations. We will use the 
intrinsic depletion zero points and slopes of the linear 
depletion relations to summarize each galaxy’s depletion 
relation. While we do not have a quantitative ” depletion- 
relation-relation” to describe variations between ensem¬ 
bles of pencil beam observations of different galaxies, we 
can use the changes we see from the MW to the LMC 
to the SMC as a qualitative depletion-relation sequence, 
at least to metallicities above that of the SMC. Starting 
from the MW and moving towards lower metallicity, in¬ 
trinsic depletion zero points decrease in magnitude, while 
depletion slopes remain approximately unchanged. The 
intrinsic depletion zero points of P and Zn, whose mag¬ 
nitude is already consistent with zero, do not change. 

Comparing depletion relations between galaxies takes 
a (generally) prohibitive amount of information. Com¬ 
puting an abundance requires absorption-derived hydro¬ 
gen columns, computing an intrinsic depletion requires 
galaxy-specific reference abundances (see Sec. 3. 3), and 
computing a depletion relation requires a representative 
sample of intrinsic depletions. In a more typical case, 
one instead has to work with a metallicity instead of a 
reference abundance for each element and less than ten 
sightlines per galaxy. This was true of depletion studies 
of the MCs, remains true of depletion studies of other 
local galaxies (e.g. Welsh & Lallement 2013; James et al. 
2014), and is, in a sense, intrinsically true of DLA and 
sub-DLA studies. 

A common exploratory technique in depletion studies 
of small samples involves making a comparison between 
a new measurement and some typical MW values for the 
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Fig. 15. — A comparison between abundances relative to solar in 
DLAs (data points, Vladilo et al. 2011, Rafelski et al. 2012) and 
the diffuse neutral media of the MW, LMC, and SMC (solid black, 
red, and blue lines, MW values from Jenkins 2009, LMC and SMC 
values from Section 3.4 of this work). 

abundance of a volatile element and the abundance ra¬ 
tio between the same volatile element and a refractory 
element. If one assumes that the MW depletion relation 
applies in the new measurement’s galaxy, one can esti¬ 
mate each element’s depletion and reference abundance. 
This type of exploratory analysis, if applied to the LMC 
or SMC, would be incorrect. We show an example of this 
type of exploratory analysis in the top panel of Figure 
14. The lines in this panel are the MW depletion rela¬ 
tion shifted to the metallicity of the MW (black), LMC 
(red), and SMC (blue). The diagonal portion of each line 
corresponds to a depletion strength (F*) range of 0 (top 
right) to 1 (bottom left), and includes what Savage & 
Sembach (1996) refer to as ’’warm disk” and ’’cool disk” 
depletion patterns. The vertical line extending upwards 
from F* = 0 corresponds to the ’’diffuse halo” depletion 
pattern range from the same work. These metallicity- 
shifted MW depletion relations are inconsistent with our 
observations. 

In the bottom panel of Figure 14, the lines are our MC 
depletion relations from 3.4. Examining these lines gives 
us a qualitative visual version of the earlier statement 
that intrinsic depletion zero points become smaller with 
decreasing metallicity. On a volatile abundance versus 
refractory-to-volatile ratio plot, a galaxy’s intrinsic de¬ 
pletion zero point corresponds to the top right end of 
that galaxy’s representative line. Because, at least in 
these three galaxies, the intrinsic depletion zero point be¬ 
comes less negative with decreasing metallicity, the top 
right end of each galaxy’s line is shifted up and to the left 
as the metallicity decreases from solar (the MW) to 1/2 
solar (the LMC) to 1/5 solar (the SMC). The fraction of 
metals in dust decreases with decreasing metallicity. 

5.4. The Magellanic Clouds and damped Lyman-a 
systems 

DLAs are reservoirs of neutral gas with H I column 
densities greater than 2 x 10^° cm“^. For a review of 
DLA properties, see Wolfe et al. (2005). Metal abun¬ 
dances in DLAs range from approximately solar to less 


than 1/100 of solar, making the DLA population a valu¬ 
able resource for understanding cosmic chemical evolu¬ 
tion. The conversion between observed and intrinsic 
abundances is a long standing problem in DLA stud¬ 
ies, since one observes a single sightline and needs to 
deconvolve the effects of metallicity, abundance varia¬ 
tions (especially a-enrichment), and depletion (Pettini 
et al. 1990; Sembach et al. 1995; Kulkarni et al. 1997; 
Prochaska & Wolfe 2002; Meiring et al. 2006; Rafelski 
et al. 2012; Vladilo et al. 2011; De Cia et al. 2013). 

The depletion part of this deconvolution problem 
is usually approached using some variation of the 
metallicity-scaled MW technique described in the pre¬ 
vious section (see, e.g. Rafelski et al. 2012; De Cia et al. 
2013). Now that we have depletion relations from multi¬ 
ple galaxies over a range of metallicities, we are no longer 
limited to just the MW. While we do not have a quan¬ 
titative model for the evolution of galaxy-wide depletion 
relations, we do have the previous section’s qualitative 
picture. Figure 15 compares the measured Fe-to-Zn and 
Si-to-Zn ratios and Zn abundances in DLAs with the 
corresponding ratios and abundances derived from the 
depletion relations of the MW, LMC and SMC. We first 
examine DLAs that have more positive Zn abundances 
or more negative refractory-to-volatile ratios than the 
SMC. The evolution of the Si-to-Zn ratio with the Zn 
abundance in our galaxy sample and in this DLA sub¬ 
sample appears to be the same. The upper envelope of 
the Fe-to-Zn ratio is perhaps slightly less negative in this 
subsample of DLAs than in our galaxies, but the separa¬ 
tion of the upper and lower envelopes along the direction 
of intra-galaxy dust evolution (i.e. parallel to the galaxy 
lines) is the same as the extent of the galaxy lines. If we 
interpret a single DLA to be a pencil-beam-sized region 
in a larger system that can contain a variety of neutral 
medium conditions, the chemical evolution of the neutral 
ISM in these systems appears to be quite similar to that 
of the neutral ISM in the MCs or MW. 

Interpreting the measurements with sub-SMC Zn 
abundances is more difficult. Applying a metallicity- 
scaled version of the MW depletion relation does not give 
us an adequate description of gas-phase abundances in 
the MCs, so metallicity-scaling the SMC depletion rela¬ 
tion may not necessarily give us an adequate description 
of gas-phase abundances in low-metallicity DLAs. The 
shape of the upper and lower envelopes of the DLA abun¬ 
dances and abundance ratios in Figure 15 broadly sug¬ 
gests that the concurrent metallicity and dust-to-metals 
ratio evolution that applies from the MW to the SMC 
continues to apply to even lower metallicities. By ap¬ 
plying the previous paragraph’s interpretation of a DLA 
as a single sightline through some system and assuming 
that the metallicity and dust-to-metals ratio of these sys¬ 
tems evolves in a consistent way, one could come up with 
a more quantitative version of this statement. We leave 
this exercise for the future. 

5.5. Gas-to-dust ratios in the Magellanic Clouds 

In Section 4.2, we computed GDRs for the DNM of 
the Magellanic Clouds from depletions. These GDRs are 
summarized in Figure 10. The DNM of each MC contains 
(spatially separate) gas with a range of overall depletion 
levels. In the figure, this level increases from left to right, 
as the amount of dust per amount of gas increases. The 
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TABLE 6 

GDR MEASUREMENTS IN THE LMC AND SMC 


Location 

GDR 

Reference 

SMC DNM 

1600-2400 

Tumlinson et al. (2002) 

SMC Wing 

480-720 

Gordon et al. (2003) 

SMC Bar 

850-1275 

Gordon et al. (2003) 

SMC DNM 

3000-4500 

Bot et al. (2004) 

SMC, global average 

700 

Leroy et al. (2007) 

SMC Tail 

850-1550 

Gordon et al. (2009) 

SMC, CO-detected regions 

150-600 

Leroy et al. (2011) 

SMC, global average 

1300-2180 

Gordon et al. (2014b) 

SMC DNM 

480-2100 

This work 

LMC DNM 

400-600 

Tumlinson et al. (2002) 

LMC DNM 

200-300 

Gordon et al. (2003) 

LMC supershells 

450-675 

Gordon et al. (2003) 

LMC, global average 

270-405 

Bernard et al. (2008) 

LMC, CO-detected regions 

150-250 

Leroy et al. (2011) 

LMC, global average 

300-400 

Gordon et al. (2014b) 

LMC DNM 

190-565 

This work 


Note. — Some of the sources give GDRs as factors of the MW 
GDR. These relative GDRs were converted to absolute values by 
assuming an MW GDR of 100-150. 


GDRs of the least depleted sightlines in the LMC and 
SMC DNM are 455-565 and 1540-2065; the GDRs of the 
most depleted sightlines are 190-245 and 480-595. We 
have compared these GDRs with values from the litera¬ 
ture and found no significant discrepancies. 

The literature values are given in Table 6. All of 
their uncertainties have some overlap with the our GDR 
ranges. This does not mean that our GDR for any sin¬ 
gle depletion level is accurate. All of the values in the 
table other than ours are averages over different DNM 
conditions or, in some cases, multiple ISM phases. An 
average over a (possibly) wide range of GDRs with un¬ 
known proportions of mass at each GDR does not contain 
enough information to draw that sort of conclusion. The 
decrease in GDR from averages over more diffuse regions 
to averages over denser regions is heartening, though not 
conclusive. 


can be accurately described using a linear depletion re¬ 
lation of the type defined for the MW by Jenkins (2009), 
but with different linear coefficients. 

We find that the ISM properties of the LMC are, up to 
metallicity, quite similar to those of the MW while the 
ISM properties of the SMC have some significant differ¬ 
ences. When we compare gas-phase abundances in the 
LMC and SMC with gas-phase abundances in damped 
Lyman-a systems, we find a considerable amount of over¬ 
lap. Within this overlap region, the gas-phase abundance 
ranges of the MCs and DLAs are very similar. Regard¬ 
less of whether this similarity is meaningful, we caution 
against assuming that MW depletion patterns apply to 
low-metallicity systems. 

Below, we list a number of more specific results. 

1. The minimum Si depletion in the SMC is consistent 
with zero while the Fe depletion along the same 
sightline is substantial, suggesting that the compo¬ 
sition of some dust in the MW and SMC differs. 

2. The volatile elements P and Zn have non-zero de¬ 
pletions along most of our sightlines through both 
MCs. The maximum P and Zn depletions in the 
MW, LMC, and SMC are the same to within the 
uncertainties. 

3. Si, Cr, and Fe are systematically less depleted in 
the SMC than in the LMC or MW. 

4. We find GDR ranges of 190-565 in the LMC and 
480-2100 in the SMC. These ranges are broadly 
consistent with GDR values from the literature. 

5. The GDRs of pencil-beam-sized parcels of neutral 
medium vary by a factor of 2 in the MW and LMC 
and by a factor of 5 in the SMC. The GDR varia¬ 
tions within each galaxy suggest some combination 
of dust evolution in the neutral medium and rapid 
cycling of interstellar matter between the neutral 
medium and dense molecular clouds. 


6. SUMMARY 

In order to study the gas and dust content of the 
LMC and SMC, we have obtained and analyzed 16 HST 
COS NUV spectra and analyzed 85 archival FUSE FUV 
spectra towards the Large and Small Magellanic Clouds. 
From these spectra, we have measured P II, Zn II, Si II, 
Cr II, and Fe II column densities. We have combined 
these measurements with H I and H 2 column densities 
and photospheric abundances of LMC and SMC stars 
from the literature to compute intrinsic elemental de¬ 
pletions, solid-phase elemental abundances, and GDRs. 
We see large variations in gas-phase abundances, deple¬ 
tions, solid-phase abundances, and GDRs from sightline 
to sightline within each galaxy. The depletion variations 
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GO-13004 was provided by NASA through grants from 
the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astron¬ 
omy, Inc., under NASA contract NAS 5-26555. 


APPENDIX 

ANALYZING A SPECTRUM 

The purpose of our model of a spectrum is to measure the total column density of each species in the 

observation-appropriate Magellanic Cloud. We marginalize over, meaning integrate over the posterior probability 
distribution of, every other model parameter. The other important high-level quantities are the (ISM absorption-free) 
continuum flux as a function of wavelength c(A) and the distribution of as a function of Doppler velocity A^'(u). 
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To convert these quantities to observables, we first use the N^^{v) to compute a transmittance function T(A). Then, 
we convolve the product of the continuum and transmittance with an instrumental line-spread function (LSF) T(A, A') 
to get a model flux /(A). This model flux is the mean of a multivariate normal (MVN) distribution, from which we 
assume the data y{X) have been drawn. The covariance matrix E of this MVN distribution is the sum of a diagonal 
measurement uncertainty matrix Emeas and a not necessarily diagonal continuum uncertainty matrix Econt- 

We fit for the ISM absorption-free continuum emission using Gaussian process (GP) regression and prediction 
(Rasmussen & Williams 2005). The regression and prediction were done using the GP code george (Ambikasaran 
et al. 2014). Our main motivation for using a GP rather than, for instance, splines is that the posterior predictive 
distribution of a GP conditioned on observations can be used to exactly compute the predicted values’ covariant 
uncertainties. This covariant continuum uncertainty is an important part of our model’s error budget, particularly for 
high SNR measurements. 

We run the GP regression on a subset of the data that we find, by visual inspection, to be free of ISM absorption, 
then compute the predicted mean and covariance matrix of the full dataset. This subset comes as close to the ISM 
absorption as possible in order to minimize the size of the prediction region. For the GP kernel, we use the sum of a 
square exponential kernel, a quadratic kernel, and the measurement uncertainties of the ISM-free subset. We fix the 
lengthscale of the square exponential kernel to be 6 A and fit for the rest of the kernel parameters. 

The dashed line in each panel of Figure 4 is an example of a continuum fit. The uncertainty of the full spectrum 
model, which includes the continuum uncertainty, is shown as a gray region in this figure. 

Instead of explicitly modeling the transmittance T(A), we instead model the species’ column density distributions 
as a function of Doppler velocity and use a deterministic non-linear transformation to convert that to T{X). 
The column density distributions are split into spatio-kinematic groups Sj , which are velocity ranges that we a priori 
assign to the MW, an I/HVG, or a MG, and the Sj are further split into components t = 1,... ,Ts ■ Each component 

■ . s ■ . .X 

has a central velocity \ a, width and, for each species X^, a fraction Ug. of the group-level column density of 
X^. The column density of species X^ in component t is then x note that the sum of the a for a single 

species over all of the components of a group are equal to 1. 

The column density distribution of a species is the sum of the group level distributions, which are in turn the sum 
of the component level distributions. The column density distribution of a single component is 






exp 



and the full column density distribution of a single element is 




3 = 1 t=l 


(Al) 


(A2) 


We use this two-level formalism because we believe that the groups are physically meaningful while the components 
may or may not be. The groups are physically meaningful because they correspond to different galaxies with reasonably 
well-separated systemic velocities. We believe that column density that has been assigned to a given group most likely 
actually arises in the corresponding galaxy. The components are not necessarily meaningful because, firstly, our 
instrumental LSFs are broader than the expected “individual” components, secondly, components that are distinct 
in velocity are not necessarily distinct spatially, and thirdly, the fragmentation of a continuous density into discrete 
subcomponents is almost always not unique. One consequence of the third point, particularly when the first point 
applies, is that we cannot even define a unique number of components. This is especially true when our observations are 
a product of the transmission function, which is a non-linear transformation and superposition of the column density 
distributions, with the uncertain continuum. Instead of trying to choose a fixed number of unphysical components, we 
marginalize over it. This means that the dimension of our parameter space is not constant. 

We generate samples from the posterior probability over this variable-dimensional parameter space, using the 
Reversible-Jump Markov Chain Monte Carlo (RJMCMC) algorithm (Green 1995). RJMCMC is a type of MCMC that 
can be applied to variable-dimensional parameter spaces. It works by considering each fixed-dimensional parameter 
space as a subspace of some larger, overarching parameter space. One then defines different Markovian steps for moving 
around within each fixed-dimensional subspace and between ‘adjacent’ fixed-dimensional subspaces. Further details 
of this framework are beyond the scope of this appendix and can be found in (Green 1995). We use these samples to 
build posterior probability distributions for each species’ column density in the LMC or SMC and the model flux at 
each wavelength. The 16-th and 84-th percentiles of the column density posterior probability distributions are shown 
as error bars in Figure 5. The 5-th, 16-th, 84-th, and 95-th percentiles of the model flux marginals are shown as shaded 
regions in Figure 4. These marginals include measurement uncertainty, continuum uncertainty, and the possibility of 
observationally similar but physically different velocity component structures. 
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CLOUDY PARAMETERS 

Here, we provide the complete list of input parameters needed to reproduce the photionization calculations described 
in section 3.2. These calculations were done using version 13.02 of the Cloudy photoionization code (Ferland et al. 
2013). All of the models are of a constant-density cloud illuminated by a star 10^® cm from the cloud surface. We 
compute models with cloud volume densities of 0.1, 1, 10, and 100 cm“^. The cloud composition is given by the 
‘abundances ism’ command, scaled to the appropriate metallicity (0.5 for the LMC, 0.2 for the SMC) using the 
‘metals grains’ command. We run each cloud model for a BO-type and an 02-type illuminating star of the appropriate 
metallicity using stellar atmosphere models from Lanz & Hubeny (2003). 

In addition to the illuminating star, the incident radiation field includes the cosmic microwave background and the 
local MW radiation field (‘table ism’) scaled by a factor of 4 in the LMC (Bernard et al. 2008) and 10 in the SMC 
(Sandstrom et al. 2010). We include the local cosmic ray background using the ‘cosmic rays background’ command. 

ELEMENTAL ABUNDANCE META-ANALYSIS 

In this appendix, we describe a multilevel model for the mean elemental abundance /r of an element from a collection 
of possibly biased studies with some repeated observations. A description of multilevel models and their application 
to meta-analysis problems can be found in Chapter 5 of Gelman et al. (2013). We can split the model into four blocks 
- population level, study level, individual star level, and study-star level. 

The population level includes the mean abundance /r, the intrinsic star-to-star variance Ustan and the between-study 
variance a study- 


y - 

Normal (/i = -|- Zmc, = 0.5^) 

(Cl) 

^star 

Exponential (A = 5) 

(C2) 

^study 

Exponential (A = 5), 

(C3) 


where p,© is the solar abundance of the relevant element, Zmc is the metallicity of the relevant Magellanic Cloud, and 
the notation “a; ~ Distribution”, as in /r ~ Normal, means that the specified distribution is our prior for x. The study 
and star levels include the bias terms Agtudyj, for study-wide effects, and AstarkJ for star-to-star variations: 

Astudy, ^ Normal (/r = 0, = a'^study) 

‘^star, ^ Normal {fd = 0, 

The study-star level consists of observations yj^k ^^d observational precisions ^ 

'^jik Gamma = C^studyj ? ^ — ^studyj^ 

yj^k ~ Normal (y = y, + Astudy^ + Astark,^-^ = 1 

^studyj — I I ? Pstudyj = 2 ’ 

V J ^Tstudyj 

where y-Tstudy- is the typical precision of the measurements in study j and cTr^tudy is the dispersion in the typical 
precision. 

We implemented this model using PyMC (Patil et al. 2010) and generated samples from the posterior probability 
distribution using the No-U Turn Sampler. We ran each element’s MCMC chain for 500 burn-in steps and 5000 kept 
steps starting from the maximum a posteriori value, which, by visual inspection, appeared to be a long-enough burn-in 
phase for convergence. 


(C4) 

(C5) 

(C6) 

da) (C7) 

(C8) 
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